Additive noise quenches delay-induced oscillations 



Jeremie Lefebvre 1 '^] and Axel Hutt 2 

1 Departement de Neurosciences Fondamentales, University of Geneva, 
CMU , 1 Rue Michel Servet, 1211 Geneva 4, Switzerland 
B INRIA CR Nancy - Grand Est, 615 Rue du Jardin Botanique, 54602 Villers-les-Nancy Cedex, France 

(Dated: March 15, 2013) 

Noise has significant impact on nonlinear phenomena. Here we demonstrate that, in opposition 
to previous assumptions, additive noise interfere with the linear stability of scalar nonlinear sys- 
tems when these are subject to time delay. We show this by performing a recently designed time- 
dependent delayed center manifold (DCM) reduction around an Hopf bifurcation in a model of 
nonlinear negative feedback. Using this, we show that noise intensity must be considered as a bifur- 
cation parameter and thus shifts the threshold at which emerge delay-induced rhythmic solutions. 
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INTRODUCTION 

Noise has a strong influence on the dynamics of systems subject to time-delayed retroaction. Random fluctuations 
impact laser dynamics 1J, gene regulatory networks 2 , postural and mechanical control [3], which all exhibit 
retarded self-interactions. In neuroscience, the propagation and time evolution of neural activity has been shown 
to vary according to combined action of sustained stochastic driving and recurrent feedback, key components 
of brain circuits @]. As such, proper characterization of systems away from the stationnary regime is essential 
in order to understand how these elements conspire and give rise to the complex phenomena we observe in experiments. 

Systems evolving under mixtures of noise and nonlinear feedback express divergent properties from their de- 
terministic counterparts, where noise has been shown to provoke or prevent bifurcations 5j. However, such 
noise-induced stability shifts in scalar systems have are thought to be caused only by a subgroup of stochastic 
influences, namely multiplicative or parametric noise. This specific type of noise has indeed been shown to be 
the source of state-transitions in nonlinear systems, while its implication on the stability of delay equations is a 
well-documented phenomenon [5HH] • 

In the past numerous studies have determined or assumed that additive noise, that is stochastic fluctuations 
independent of the state variable, has no impact on the stability of nonlinear systems [3], even in the presence of 
delays |10j . This implies that the deterministic expression of the dynamics fully exposes the mean behavior of an 
otherwise stochastic system. This view has since been challenged in mean-field type systems [111 112] and physical 
pattern forming systems [13] . More recently, it has been demonstrated that additive noise stabilizes delayed nonlinear 
systems near non-oscillatory instabilities [IU[T5]. In this Letter, we show that additive stochastic fluctuations cause 
a shift in the eigenspectrum in a nonlinear scalar delayed differential equation(DDE) in the vicinity of an Hopf 
bifurcation and de facto shapes delay-induced rhythms. By doing so, we state that despite the current belief, additive 
noise changes significantly the stability of delayed oscillatory nonlinear systems and that the additive noise intensity 
has to be considered as a bifurcation parameter. 

To demonstrate this, we adopt the following strategy. We first derive a time-dependent delayed center mani- 
fold (DCM) representation of a generic scalar DDE unfolded around a supercritical Hopf bifurcation. Then, adiabatic 
elimination is performed based on the characteristic separation of time scales involved in the vicinity of a non- 
hyperbolic fixed point. We show that the explicit time-dependence of the DCM reveals noise-dependent components 
reponsible for the depletion of rhythmic solutions for negative nonlinearities, and consequently exposes the impact 
of additive noise on the system's equilibrium and its stability. We then apply our theory to reconstruct the associ- 
ated stochastic bifurcation diagram of the ensemble averaged noisy solutions, and validate it via numerical simulations. 
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MODEL 

The work is motivated by non-homogeneous retarded equations with the generic form 

X(t) = -X(t)+ I f(x(t-T)) + I(t), (1) 

exhibiting a smooth nonlinearity f(x) and some additive random fluctuation I(t), in the ItA' sense. Let us consider 
the expansion of Eq. ([!]) around the stationary state x — jf(x ) up to third order, which yields the polynomial 
system 

u(t) =-u + r]U T + nu 2 T + vu\ + V2D£(t) , (2) 

where r\ — jf'(x ), k — jf"(x )/2\ and v = 7/"' (x Q )/3! are non-zero constants, with the small deviations from the 
stationary state u — u(t) and u T — u(t — r). Here, the additive random uncorrelated fluctuations are defined by 
I(t) — \/2Z)£(i) with the noise intensity D > such that (£{t')£(t)) = 8(t — t') where (•) denotes the ensemble average. 

Figure [I] shows significant differences between the dynamics of Eq.([2]) for D = and D > close to the bi- 
furcation point. While virtually masked by local fluctuations in single trial trajectories, the investigation of the 
ensemble average reveals that additive noise stabilizes the oscillations by preventing the bifurcation. This is caused 
by a shift of the eigenspectrum performed by noise and delay. The impact of noise decreases as |e| increases i.e. 
away from the bifurcation, the time-scale separation vanishes and the system's ensemble average becomes somewhat 
passive to addive noise. 
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FIG. 1: Impact of additive noise on the stability of delay induced oscillations. From panel a) to c), the system Q 
undergoes a supercritical Hopf bifurcation as the control parameter e = rj — r\ c is varied from 0.25 to —0.25 and crosses the Hopf 
threshold located at e = 0. Slightly above the bifurcation in panel b), one observes deviations between the noise-free solution 
u(t)(red curve; D = 0) and the ensemble averaged solution (u(t)) (dark grey curve) over independent realizations of the noise. 
While the noise-free system exhibits stable oscillations of fixed amplitudes, noise prevents the bifurcation and mean solution's 
amplitude decreases (panel b)). The deviation between noisy and noise-free dynamics vanishes as one moves away below (panel 
a)) or above (panel c)) the bifurcation threshold. This effect is not easily observable amongst the single-trial realizations 
(light gray). DCM corrected solutions capture the effect of the additive noise on stability and reproduces the dynamics of 
the ensemble average (dashed curve, almost undistinguishable from the ensemble average). Parameters are D = 0.00001, 
7 = —0.05, t = 12.0, Tj c — —1.03 and f(x) = (1 + exp[— 60 • x])~ . The time integration followed a Euler-Maruyama scheme 
with 10 4 iterations and time step dt — 0.1. 



UNFOLDING AROUND THE HOPF BIFURCATION 

In order to explain this novel phenomenon, the subsequent paragraphs at first neglect noise and follow the standard 
analysis of DCM reduction for deterministic delay equations 16 -19 . Analyzing the stationary state's linear stability, 
we consider just the linear terms in Eq. (|2| 



u(t) = —u + rju T 



(3) 
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and focus on the Hopf instability assuming u(t) = uooxp(At), A £ C, uq =const. Then there is a critical A c = ±iw c 
that satisfies 

X c = -1 + Vc e x ^ , (4) 

where w c = yfrfc — 1 is the corresponding Hopf frequency. Introducing the control parameter e = r\ — r/ c , it is 
convenient to re-write Eq. ([2| as 

u(t) = L(u, 77 C ) + F(e, u) + V2D£(t), e(t) = , (5) 
with L(u, ?7 C ) = — m + ?; c u r and F(e, u) = eu T + ku^. + vu^.. In this formulation, the Hopf bifurcation occurs for e = 0. 



DCM REDUCTION NEGLECTING NOISE 



Equation (|3| and the characteristic equation Q define eigenvectors = {4>\{9), 4>2{&)) where 9 is a parametrization 
of the delay i.e. — r < 9 < 0. They span the center subspace and are determined by ±iw c 4>j(0) = L(<pj,rj c ) — 
-M°) + Ve4>j(- T )\3 = (M) yielding 

S(0) = (0x(0), &(0)) = (e™< 8 , e"™^) . (6) 

The associated adjoint basis of the center subspace ^(s) = (V>i(s), ^(s)) 1 = (de~ iWaS , de ZWcS ) J , d E C, obeys ($, \E') = 
I, where {a(0), b(8)) is a bilinear form defined by [Lunel] 

(a, b) =o(0)6(0)-/° [ a(Z-6)b(OdZ[da(0)]. (7) 



Here, [da(0)] = (— <5(#) +-q c S(9 + r))d9. According to the center manifold theorem, phase space can be split into center 
(C) and stable (S) parts i.e. C © 5, such that the stables modes converge to the center manifold. We thus have 

ut(0) = §{6)z + h{z,6) e C 

for some stable manifold h, where z = [z\(t), z% (t)] is the amplitude on the center manifold. The dynamics of Eq. ^ 
projected onto the linear center subspace is consequenly given by the order parameter equation 

z(t) = Bz + *(0)F[e, $z + h], e(t) = , (8) 

where the diagonal matrix B is defined by 

— lW r 



STOCHASTIC DCM REDUCTION 



Now we extend this standard DCM analysis close to a Hopf instability by introducing additive noise i.e. D > 0. We 
treat the additive noise similar to the nonlinearity which renders the stable and center manifolds time-dependent. 
Hence Eq. ^ reduces to 

z(t) = Bz + %(Q)F[e,<f>(d)z + h(z,9,t)} 

W2D*(0)e(i), e(t) = . (9) 

The separable ansatz for the time dependent DCM close to a non-oscillatory instability [HI [22] reads 



h(z,0,t) = h(z,9) + h t (9,t) 



(10) 
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Here h t (0,t) = (1 — P c )H(t + 9) is independent of z and evolves on the linear stable subspace. The operator P c 
projects onto the center subspace spanned by ( I > (#)[15] and H(t) £ C obeys 

H'(t) = L(H,r lc ) + V2DZ(t). (11) 



The separation ansatz in Eq. ( 10 ) provides a good approximitive description of delayed noisy dynamics near 
instabilities in regimes where $(0)z + h(z) evolves slowly with small amplitude. In the vivinity of an oscillatory 
instability this corresponds to regimes where the emerging oscillation has a small frequency and small ampli- 
tude, which occurs when r is taken large enough. This is the regime we investigate here. This choice is further in 
line with past numerical findings where the impact of noise was shown to be inversely proportional to the delay size 20] . 

The state probability density p(z,ht) of the system evolving according to Eq. Q considers h t as a time-dependent 
random variable and obeys the Fokker-Planck equation, 



dt YLK ' \ d 2 Zl 
-g^ {iw cZl + F[g(z) + h t ])p(z, h t ) | + V> 2 (0) 

D 9p2 ^ ht) - A (-iw c z 2 + F[g{z) + h t ])p(z, h) 
a 2 z 2 oz 2 



(12) 



where g(z) = $(6>)z + h(z,9) is the time-independent contribution to the center manifold. Here the dependence in 9 
and the trivial dynamics of the unfolding parameter e have been dropped for clarity. 

The characteristic time scale separation involved near non-hyperbolic fixed points implies that the time-dependent 
functional ht(t) and the noisy input £(i) are much smaller than the center modes z and retain the time-scale 



separation of the unperturbed system. By virtue of the separable ansatz in Eq. ( 10 ) for which the temporal correction 



h t is independent of z, the adiabatic elimination |13l 12 1] of the fast time dependent terms applies by performing an 



ensemble average of p{z,h t ) over ht for all 9 in Eq. (12) using the associated density q(ht). Specifically, for weak 
noise one separates the slow and fast evolution by p(z, h t ) — Po(z, t)q(h t ) and expand the nonlinearity F with respect 
to h t around the center modes. Assuming stationarity of h t and a Gaussian profile for p(h t ) with vanishing mean 
and variance a 2 = <r 2 (D) = J Q T h 2 (9, s)p{h t {9, s))ds > for a short time window T ^ 2n/w c shorter than the time 
scale of the center dynamics 

i-t+T 



F[g(z(t)) + h t (9,s)}p(h t (6,s))dsn 

eg T (z(t)) + k {g 2 T (z{t)) + a 2 ) + v (g*{z(t)) + 3g T (z(t))a 2 ) . 

with g T (z{t)) = 4>(— r)z(i) + h(z(t),—T). The resulting adiabatic Fokker Planck equation determines the temporal 
evolution of the probability density po(z, t) whose mean value z obeys the non-autonomous order parameter equation 

z(t) = Co + Bz{t) + *(0)F(e*, ${9)z(t) + h(z{t), 9)), 

(13) 

with the effective control parameter e*(D) = e + [i(D) : c (D) — kct 2 (D) and n(D) — 3v<r 2 (D). 

A comparison between the full and reduced description reveals the mechanism by which additive noise alters the 



properties of the noise-free system. To see this, note that Eq. ( 13 ) may be viewed as the reduced dynamics of the 
non-autonomous DDE 

Z(t) = c (D)+L(Z(t),r, c )+F(e*(D),Z(t)) (14) 

This equation is our major result. By direct comparison with Eq. (J5|, two effects can be seen in this system with 
cubic nonlinearity: i) Since c a ^ 0, for D > the system has a new equilibrium repositioned in phase space and, 
ii) the noise shifts the control parameter by a factor of [i and thus alters the stability of the system. Note that Eq. 
(14 1 converges to Eq. ^ when D = 0. Summarizing, the sign and manitude of the nonlinear parameters k and v in 



5 



the expanded system fully determine the qualitative impact of the additive perturbations on stability. This finding 
demonstrates the purely nonlinear nature of phenomenon shown in Fig. [TJ 

To test our theory, we have numerically investigated the behavior of the DDE in Eq. ^ with D = and D > where 
the delayed nonlinearity is taken to be negative (7 < 0). To compare the original dynamics with the DCM corrected 
system in Eq. ( 14 ) , we have chosen a regime of large delays where the center modes oscillate at a small frequency and 
where the impact of noise on the dynamics is known to be weak [20 . To determine a 2 , we have integrated numerically 
the linear system in Eq. ( |11[ ), applied the corresponding definition of h t , and computed the stationnary variance 
a 2 w (1 — 2Re(d) cos{w c t)Y{H 2 ) , which holds true for D > small. The constant d is the normalization constant of 
the adjoint basis ^f(s). 

While set in such a regime of slow oscillations, Fig. [2] shows that additive noise shifts the Hopf bifurcation to smaller 
values and thus stabilizes the system. The corrected system ( 14 ) captures very well the dynamics altered by additive 
noise and further allows to predict the position of the noise-shifted Hopf bifurcation in parameter space. Indeed, the 
DCM reduction scheme captures both qualitatively and quantitatively the effect of additive noise on the dynamics of 
the original equation in the neighborhood of the bifurcation whenever D and |e| are chosen small. It further reveals the 
characteristic quenching action of the random fluctuation on the ensemble averaged solutions experienced in presence 
negative delayed feedback. 
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FIG. 2: Delayed stochastic Hopf bifurcation diagram. Amplitudes of the nascent oscillatory solutions as a function of 
noise intensity (D) and distance from the bifurcation threshold e. The Hopf threshold is located at e = (dashed line across 
plots for comparison), a) Generic bifucation diagram in absence of noise (D = 0). b) As noise is injected additively in the 
system, the onset of stable oscillatory solutions retreats in parameter space compared to the noise-free case depicted in panel 
a), giving rise to significantly different behavior. Th e panel shows the amplitudes of the ensemble averaged solutions (u(t)) of 
Eq. (J2J . c) Amplitudes of the corrected system (14 1 capture the additive noise effect. Deviations between the noisy diagram 
and the DCM reduction prediction (between panels b) and c)) appear as D and |e| increase pointing to the approximative 
nature of the separation ansatz (jTOD . Other parameters are taken from Fig. [I] 



One may ask whether the effect revealed, i.e. the shift of the control parameter and the oscillation amplitude, is 
prominent enough to play a significant role in system's dynamics. Figure [3] shows oscillation amplitudes subjected to 
the control parameter e for two noise levels. It reveals that the oscillation amplitudes are practically indistinguishable 
for e < —0.25 or e > 0.02, i.e. away from the original deterministic bifurcation threshold. Hence the additive noise 
effect on the ensemble average of the system activity is observable near deterministic stability threshold but disappears 
far from the bifurcation point. 



DISCUSSION 



Our results reveal that additive noise perturbs the stability of scalar delayed systems and shifts the instability 
threshold in parameter space, demonstrating that noise intensity must be considered as a bifurcation parameter. 
Using time-dependent DCM reduction, we have exploited the characteristic time-scale separation emerging near 
non-hyperbolic fixed points to perform adiabatic elimination. For the cubic system considered, our approach reveals 
the effect of additive noiseon linear and cubic terms. This perturbs the stability of the systems and shifts the onset 
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FIG. 3: Bifurcation diagram for oscillation amplitudes subjected to the control parameter. The amplitudes of 
the oscillatory solutions are computed by numerical integration of the original stochastic DDE in Eq. |2| for D = 0(red) and 
D = 0.00001 (black). We observe that the two curves are well distinguishable in the neighborhood of the bifurcation point at 
e = 0, but coalesce otherwise. Other parameters are taken from Fig. [l] 



of oscillatory solutions. 

This novel finding has several implications. First, the effect suggests a general mechanism for on-line sup- 
pression of rhythmic activity by external stochastic driving in presence of inhibitory feedback, despite purely 
deterministic model parameters. We conjecture that this mechanism might provide support to experimental 
paradigms seeking to supress undesired or pathological rhythmic regimes in recurrent systems. 

The work considers the supercritical Hopf bifurcation for a delayed nonlinearity with negative gain (7 < 0). 
Although our theory is not specific to this case, we note that in the vicinity of a subcritical instability, additive noise 
causes the solutions to diverge and that, in most practical applications, the supercritical case prevails. This case 
nonetheless suggests that noise in such systems might have the opposite effect: it might enhance oscillatory behavior. 
This remains to be shown in future work. 
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